Preprint typeset in JHEP style - HYPER VERSION 



Electroweak phase transitions in the secluded 
£7(1) -extended MSSM 



Cheng-Wei Chiang 1 ' 2 and Eibun Senaha 1,3 

1 Department of Physics and Center for Mathematics and Theoretical Physics, National 
Central University, 300 Jhongda Rd., Jhongli, Taiwan 320, R.O.C. 
2 Institute of Physics, Academia Sinica, 128 Sec. 2, Academia Rd., Nankang, Taipei, 
Taiwan 115, R.O.C. 

3 Physics Division, National Center for Theoretical Sciences, 101, Sec. 2 Kuang Fu Rd., 
Hsinchu, Taiwan 300, R.O.C. 



E-mail: chengwei@phy . ncu.edu. tw, senaha@phys.cts.nthu.edu.tw 



Abstract: The electroweak phase transition (EWPT) in the secluded-C/(l)'-extended 
MSSM (sMSSM) is studied. Using the effective potential at zero and finite temperatures, 
we search for the non-MSSM-like EWPT in which the light stop mass is larger than the 
top quark mass. Scanning the parameters relevant to the EWPT, the upper limits of the 
Higgs boson masses, which are consistent with the strong first order EWPT, are derived. 
For the lightest CP-even and -odd Higgs bosons, we find vnjj x < 160 GeV and < 250 
GeV, respectively. In the sMSSM, the tree-level CP violation is possible by the complex 
soft supersymmetry breaking masses. It is observed that such a CP-violating effect does 
not spoil the strong first order EWPT for the typical parameter sets. 
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1. Introduction 

To explain the baryon asymmetry of the Universe (BAU) is one of the challenging problems 
in particle physics and cosmology. From the latest cosmological observations, the BAU is 
found to be §. 

— = (4.7 - 6.5) x 1(T 10 (95% C.L.), (1.1) 

where hb is the difference between the number density of baryons and that of antibaryons, 
and n 7 denotes the photon number density. If there is an inflation in the early Universe, any 
primordial BAU would be washed out. Therefore, the BAU must arise dynamically after 
the inflation. In order to generate the BAU from a baryon-symmetric universe, it is required 
that [Q] (1) baryon number (B) violation, (2) C and CP violation and (3) departure from 
thermal equilibrium. The last condition is not mandatory if the CPT theorem does not 
hold. In principle the standard model (SM) can satisfy these three conditions. However, 
it turns out that the CP-violating phase in the Cabibbo-Kobayashi-Maswaka matrix || is 
way too small to generate the observed BAU 0] and a strong first order electroweak phase 
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transition (EWPT) cannot be realized with the viable Higgs mass, > 114.4 GeV Q, 
rendering condition (3) infeasible ||. So far, many baryogenesis scenarios beyond the SM 
have been proposed. Among them electroweak baryogenesis || is an attractive idea since 
it exclusively relies on electroweak physics which is testable at colliders or cosmological 
observations in the near future. The minimal supersymmetric standard model (MSSM), 
which is one of the well-motivated candidates for new physics at TeV scale, may have a 
window to a successful baryogenesis |7|, |8|, ||. However, due to the strong experimental 
constraints the viable region seems to be quite limited. 

A simple way to extend the MSSM is to introduce a gauge singlet (S) into the su- 
perpotential. Several versions of the singlet extended MSSMs have been proposed: the 
Next-to MSSM (NMSSM) © 0], the nearly MSSM (nMSSM) or the minimal non- 
MSSM(MNMSSM) @, the [/(l)'-extended MSSM (UMSSM) ||, ||,[l| and the secluded 



[/(l)'-extended MSSM (sMSSM/S-model) 0, [17|, gj, |g |2| etc. In this class of the mod- 
els, there is no fundamental fi parameter in the superpotential, and it is generated after 
S develops its vacuum expectation value (VEV), i.e., fi e g = Xvs where A denotes the di- 
mensionless coupling and vs is the VEV of S. It thus gives a solution to the so-called \x 
problem. 

It is known that the Higgs singlet which couples to the Higgs doublets can play a role 
in strengthening the first order EWPT, and a light stop is not necessarily lighter than top 
quark in contrast to the MSSM case. A variety of patterns of the EWPT are possible in the 
NMSSM ill. Detailed studies of the EWPT in the nMSSM can be found in Refs. M, M. 



For the analysis of the EWPT in the UMSSM, see, for example, Refs. |p4| . The possibility 
of electroweak baryogenesis in the sMSSM is outlined in the letter paper [18], and recently 
its full paper has come out ]2(|. It is demonstrated that the electroweak baryogenesis 
is successful in this model. Although a number of studies on the EWPT in the singlet- 
extended MSSM can be found in the literature, it is still not clear how heavy Higgs boson 
can be consistent with the strong first order EWPT in some models. Such a mass limit 
would be indispensable for a test of EW baryogenesis scenario at colliders. 

In this paper, we examine the sMSSM EWPT with/without CP violation in the wider 
parameter space which has not been probed in Refs. (l^, |2C|| . Since there are the singlet 
contributions in the Higgs potential, the stability of the Higgs potential is not manifest. In 
our analysis, we require that the prescribed EW vacuum should be a global minimum at 
zero temperature. This vacuum condition can make the allowed region quite limited [19]. 
To investigate the properties of the EWPT, we use the one-loop effective potential at zero 
and nonzero temperatures taking the contributions from the Z and W bosons, the third 
generation of quarks and squarks into account. Under the theoretical and experimental 
constraints, we exclusively search for the non-MSSM-like EWPT by scanning the relevant 
parameters and work out the upper limits of the Higgs boson masses which are consistent 
with the strong first order EWPT. 

In the sMSSM, owing to the soft supersymmetry (SUSY) breaking mass terms, it is 
possible to realize both explicit and spontaneous CP violation at the tree level. The effect 
of such a CP- violating phase on the strength of the first order EWPT and the Higgs mass 
spectrum are investigated. 



- 2 - 



In general, there are ten order parameters if spontaneous CP violation exists. It is a 
non-trivial task to investigate such an EWPT thoroughly, and the numerical calculation is 
extremely time consuming. We will focus exclusively on both the CP-conserving case and 
the explicit CP-violating case in which the number of order parameters is reduced to six. 

The paper is organized as follows. In Sec. ^, we briefly explain the model and define 
the notations. Section || is devoted to the tree-level analysis. The approximate formulae 
of the Higgs spectrum are presented. We discuss the vacuum condition, which turns out 
to be the severest theoretical constraint in this study. In Sec. ||, we argue the EWPT 
qualitatively in some detail. To get some idea about the sMSSM EWPT, we outline the 
nMSSM EWPT by showing the analytical formulae. In Sec. [Si, the numerical results are 
presented. Sec. ^ contains the conclusions and discussions. 



2. The model 

The sMSSM is one of the singlet-extended MSSM models, which can be regarded as an 
effective theories of some unification theory such as string theory. The symmetry of the 
model is SU(3)c x SU(2)i x U(l)y x U(1)'q,, where the extra U(l) is a remnant of the 
larger symmetry of the UV theory. The Higgs sector comprises two Higgs doublets (Hd, H u ) 
and four Higgs singlets (S, S±, S2, S3). Among the singlets, the so-called secluded singlets 
(Si,i = 1,2,3.) play an essential role in ameliorating the severe experimental constraints 
on the Z' boson fl(| . 

It is desirable to have U{\)' charges (Q's) chosen to make the model anomaly free. 



To this end, exotic chiral supermultiplets are generally required 14, g^, £6|. For our 
purpose, we assume that they are heavy enough not to affect the phenomenology at the 
electroweak scale. Neither will we address the gauge coupling unification issue here as it 
requires the knowledge of full particle spectrum in the model. Instead, we focus exclusively 
on the Higgs sector to discuss the EWPT. Here, the so-called Model I sMSSM is considered, 
which can accommodate CP violation at the tree level. In Model I, the U(l)' charges of 
the Higgs fields must satisfy 

Q Hd + Qh u + Qs = o, Qs = -Qsi = ~Qs 2 = \Qs s - (2.1) 

The superpotential (W) includes the following trilinear terms: 

W B -eijXSHjHl - XsS^Ss , (2.2) 

where H^ u , S,Si,i = 1,2,3, are the chiral superfields, A and Xs are the dimensionless 
couplings. Due to the U(l)' symmetry, the S n (n 6 Z) terms are forbidden. Therefore, we 
are not bothered by the domain wall problem, which can be induced by the spontaneous 
breaking of the discrete symmetries. In general, the terms of the form SS\, SS2, SfS^ and 
S$ S3 are allowed under the charge assignments (^]l|). However, we will not consider such 
quadratic terms since they may reintroduce the (i problem. Neither will we consider the 
cubic terms for simplicity. As we mentioned in the Introduction, once the Higgs singlet S 
develops its VEV, the effective fj, term is generated by fj, e g = X(S). Therefore, the scale of 
H e g is not completely arbitrary but is fixed by the soft SUSY breaking parameters. 
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The Higgs potential at the tree level is given by 



Vo = V F + V D + V soit , (2.3) 

with 

V F = |A| 2 {|e^<FJ 2 + \S\ 2 (^ d + <S>l$ u )} + |A 5 | 2 (|SiS 2 | 2 + \S 2 S 3 \ 2 + |S 3 5i| 2 ), (2.4) 
V D = ^±^{^ d - + 9 4\^ u \ 2 

o Z 

+ 9 4 (Qn d & d *d + Qh u ^u + Qs\S\ 2 + J2 Qs.\S*\ 2 ) \ (2-5) 



Koft = + m 2 <S>l$ u + m||5| 2 + J^m|j5, 



i=l 
2 



8=1 



-(eijAA A S'$^i + \sA\ s SiS 2 S 3 + m 2 SSl SSi + m 2 SS2 SS 2 + m 2 SlS2 S\S 2 + h.c), 

(2.6) 

where <?i and ^ are the SU(2), U(l) and i7(l)' gauge couplings, respectively. We will 
take g[ = y/5/3gi as motivated by the gauge unification in the simple GUTs. The soft 
SUSY breaking masses rn 2 SSi and ni 2 SS2 are introduced to break two unwanted global U(l) 
symmetries, and rng iS 15 nee ded for explicit CP violation (ECPV). 

Here, we will follow the notation given in Ref. The Higgs VEVs and their fluctu- 
ation fields are parameterized as 

= e * ( ^ {Vd \ hd + md) )■ *> = ( ' , tl + ■ > V < 2 - 7 > 

e id s e i8 Sl 

S = —^-(vs + h s + ia s ), Si = -j=-{v Si + h Si +ia Si ), i = 1 - 3, (2.8) 



where " = y^ + ^u - 246 GeV at the vacuum. The nonzero 9's can bring about spon- 
taneous CP violation (SCPV). However, not all #'s are independent. Due to the gauge 
invariance, the following four combinations are physical: 

<pi = 0s + 0si, f2 = o,s + o S2 , ^ = o s + e l + e 2 , = o Sl + + e Sa . (2.9) 



The analysis of SCPV at zero temperature can be found in Ref. IS]. To accommodate 
SCPV, the lightest Higgs boson mass should be less than about 125 GeV. Although it is 
interesting to study SCPV at finite temperature, we will not pursue this possibility in this 
paper. A comprehensive study of the SCPV at the finite temperature can be found in 
Ref. fM. 
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The tadpole conditions, which are defined by the first derivative of the Higgs potential 
with respect to the Higgs fields, are given by 



1 /dV \ 2 , gf+ffj / 2 2n p v uV S , \M ,2, 2N,5l n a n 



das/ 
dV 
da Sl 
dVp 
da S2 
/ dV 
\da Ss 



where 



(2.10) 



l/dV \ 2 52+51/2 2^ o , |A| 2 2 2 fff 



V u \dhJ ~ 1112 ~ ~^ yVd ~ Uu > ~ ~ K d s; T 

(2.11) 

1 / dV \ 1 IAI 2 a' 2 

— {7^)= m s- ( R i v ^ + R ^ + Rxvdvu)— + ^(vj + vl) + ^-Q S A = 0, 
vs \ oh s I vs 2 2 



~~(^r^/ = m2 si ~ ( RlV s + ^12^52 + R\ s vs 2 vs 3 )— 

VSi \ 0< l s 1 I V Sl 



1 / 9V \ 2 r, r, X 1 

- m S2 - {R 2 vs + Ruvs! + R\ s v Sl v S3 ) : 



v S2 \dh s J "•'» ' — ' ,^ 

— (v Sl +v S3 ) + - 



1 / dV \ 2 D VS1VS2 . |Ag|% 2 , 2 x , 9? n A n 



i=i 

,2 \ T T™/™2 



_ Re(A g A As ) _ Im(A 5 ^ As ; 



(2.12) 



%^1 + 4J + ^A = 0, (2.13) 



^(4 1 +-! 3 ) + f^ 2 A = 0, (2.14) 



(2.15) 



-(TZr)= I W = 0, (2.16) 



^i«Si + ^s 2 + I\v d v u = 0, (2.17) 

hvs - I12VS2 + h s vs 2 vs 3 = 0, (2.18) 

/2W5 + Invs! + I\ s vsivs 3 = °> ( 2 - 19 ) 

/A fl VSi«Sa = 0, (2.20) 



A = Qtf d ^ + Qif u ^ + Qsv 2 s + £ Qs^Ii, (2-21) 



Ri = Re(m SSi ), h = lm{m ss .), z = 1,2, (2.22) 

flia = Re(m| l52 ), J 12 = Im(m| l52 ), (2.23) 

Re(XA x ) lm(\A\) 

Rx = ^=-, / A = ^ 7f -, (2.24) 
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and {X) is denned such that X is evaluated at the vacuum. In our analysis, the soft 
SUSY breaking masses (mf, m 2 , m|, rn 2 s _,i = 1,2,3.) are determined via the six tadpole 
conditions (g!c|)-(pl5D. Here, all the Higgs VEVs are regarded as the input parameters 



and assumed to be nonzero. 

From the tadpole conditions with respect to the CP-odd Higgs fields, it follows that 

I x = I Xs =0, h=h2^, h = -h2 — - (2.26) 

v s v s 

Therefore, there is only one physical CP-violating phase at the tree level, and we take 
Arg(m| i5a ) = 9s 1 S2 as an input. After including the one-loop contributions, the relations 
( 2.26|) will be modified, and will be discussed in subsection |5.3j . 



3. Higgs mass spectrum and vacuum conditions 

In this paper, we consider the one-loop corrections from the Z and W bosons, the third 
generation of quarks (t,b), and squarks (£1,2,^1,2)- The one-loop effective potential at zero 
temperature takes the form [27| 

ViC^.S.SiAs) = 5>|J (ing - I) , (3.1) 

which is regularized in the DR scheme, rfiA is a field dependent mass, and M is the 
renormalization scale determined by the condition (Vi) = 0. The statistical factor of each 
particle is given respectively by cz = 3, cyy = 6, q = c& = —ANc, and 2 = eg = 2iVc, 
where A r c* is the color factor. 

In principle, the Z boson can mix with the Z' boson, and their mass matrix becomes 
2-by-2. However, the mass of Z 1 boson and the magnitude of the mixing angle are strongly 
constrained by experiments. A recent analysis of the constraints on the Z 1 boson mass and 
the mixing angle can be found in Ref. |28[| . As the mixing angle is constrained to be less 
than 10~ 3 pSfl , we simply consider the case of no Z-Z' mixing, i.e., tan/3 = \jQn d IQu u - 



In addition, since the tan/3 ~ C(l) is generic in the sMSSM [16, 17, 19|, we will present 



only the case of tan [3 = 1. In this case, the secluded sector does not contribute to Eq. (|Q|). 
3.1 Higgs mass spectrum 

Due to the additional contributions coming from the P-term and P-term of U(l)', the 
mass bound on the lightest Higgs boson is significantly relaxed compared to the MSSM. 
At the tree level, it is found that 



I \|2 

m 2 Hl < ml cos 2 2/3 + l -^-v 2 sin 2 2/3 + g?v 2 (Q Hd cos 2 + Q Hu sin 2 /3) 2 



= (0 GeV) 2 + (139 GeV) 2 + (111 GeV) 2 ~ (178 GeV) 2 , (3.2) 

where we have taken tan/3 = 1, A = 0.8, Qh a = Qh u = 1 in the second line. Thus, large 
radiative corrections are not necessarily required for avoiding the LEP exclusion mass 
limits. However, we should note that because of the mixing terms between the doublets 
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and singlets, rn? H can become smaller and even negative. For the above parameter set, we 
get 



2 1 
m Hl = - 



m 2 + |A| V + 6g' 2 v 2 s - Mm 2 + 2g' 2 (3v 2 - v*)} + 4v*{r x - (|A|2 - 2g' 2 )v s } 



(3.3) 



where the mixing terms coming from the secluded sector are neglected, and CP is assumed 
to be conserved for simplicity. The parameter m 2 s appearing in Eq. fl3.3j ) is given by the 
tadpole condition for hs, i-e., Eq. (2.12). Therefore, rnjj can become negative in the large 
R\ limit. It should be noted that R\ can be re-expressed in terms of m H ±. At the tree 
level, 

2 1 / d 2 V \ 2 2R X \X\ 2 2 

m H± = —-T—^{Trr+^)=rn w + -—v s -—v. (3.4) 



sin/3cos/3\d0+<9<fc/ sin 2/3 2 

Thus, rriH 1 can be unphysical in the large m H ± for a moderate value of vs- The one-loop 



formula of m 2 H± is explicitly given in Ref. [19|. As is done there, we take m H ± as an input 
in place of 

To obtain more precise values in the Higgs mass spectrum, it is necessary to incorpo- 
rate the mixing terms between ($^, and (S\, S2, S3) sectors, and also the one-loop 
contributions (denoted by Am|). It is well-known that the dominant term of £\m 2 H comes 
from the top/stop loops: 

o Nr 8 a rnr mr 
Am H = 17T2^ m t^^^- (3-5) 

The explicit formulae of the mass matrix at the tree level are presented in Refs. ]l6|, [l9|] . 
For the one-loop expression, which is the same as the that of the NMSSM, see, for example, 



Ref. [11]. After the SU{2) Nambu-Goldstone bosons are rotated away, the mass matrix of 
the neutral Higgs bosons can be reduced to an 11-by-ll form. We diagonalize it numerically 
to obtain , where the subscript i is labeled in the ascending order of mass. 

Here, we comment on the possible range of m H ± . It is known that m H ± is constrained 
from above by the vacuum condition. Namely, if we require that the energy level of the 



electroweak vacuum v = y v d + v\ ~ 246 GeV should be lower than the origin at which 
the energy level is normalized to zero, cannot exceed some critical value. For the 



typical parameter sets, the maximal value of is found to be 0(1 — 10) TeV [19]. 

Actually, we can get a stronger constraint since the symmetric vacuum where v & = v u = 
is not necessarily located at the origin. We will discuss the vacuum condition in the next 
subsection. 

3.2 Energy levels of the vacua 

Before moving on to the analysis of the EWPT, we consider the structure of the zero- 
temperature effective potential in some detail. Although we restrict ourselves to the tree- 
level analysis for simplicity, it suffices to know the qualitative features of the Higgs potential 
in the sMSSM. 
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The effective potential at the tree level takes the form 

Vo = ^rn\v 2 d + -m\v\ + -m 2 s v 2 s + ^ ^rn 2 s .v 2 s . - Riv s v Sl - R2V S vs 2 - Ri2V Sl vs 2 

i 

D D i ^2 + 9l i 2 2\2 , l^| 2 / 2 2, 2 2, 2 2\ 

-R x v d v u v s - R\ s v Sl v S2 v Sz + 32 forf - Wj + — + + v u v s ) 

i I^S 1 ) 2 / 2 2 , 2 2 , 2 2 \ , 5? a 2 /o c^i 

+^p( u Si U 5 2 + % 2 % 3 + v S 3 V Sl ) + Y A ■ ( 3 ' 6 ) 

The EW vacuum defined by the tadpole conditions (in what follows we refer it as the 
prescribed EW vacuum) is not always the global minimum. In order to see this, we compare 
the energy levels of the EW vacuum and the symmetric vacuum. The energy levels of the 
two vacua are given respectively by 



1 1 Q 2 + Q 2 

(^o)vac = -RWdVuVS + -R\ s V Sl V S2 Vs 3 ~ 2 — 1 {vj - <) 



2 „,2\2 



I^P/„,2„,2 i „,2„,2 i „,2„.2^ l^s| 2 ^,2 ,,2 , „,2 „,2 , „,2 „,2 \ 9l 
4 



-^{v 2 vl + vjv 2 s + v 2 y s ) - ^(v 2 Si v 2 S2 + v 2 2 v% + v 2 S3 v 2 Sl ) - ^A 2 , 

(3.7) 

(V^U = (V (v d ,u = 0)) vac 



1 I \ 1 2 „'2 

T> . -._ „-_ „-_ \ As \ t-,2 „-2 i -2 -2 i -2 -2 A 1/1 

(sym) 



~Rx s v Sl v S2 vs 3 - ^(v 2 Sl v 2 2 + v 2 2 v 2 3 + v 2 3 v 2 Si ) - ^A 2 . (3.8) 



where v'a are determined by the tadpole conditions using the potential Vq , and A is 
given by A(v = v). The energy level of the symmetric vacuum can in principle become 
higher than the origin in the limit of large positive R\ s . 

The difference between the energy levels of the two vacua is 

A(y )vac = (^) (Sym) )vac " (Vo)vac 

= ~^R\VdVuVs + ^Rxsivs^SiVs-s ~ vs 1 vs 2 vs 3 ) 

- -If + + + vlvl) 

|As| 2 



4 



-2 -2,-2 -2,-2 -2 2 2 2 2 2 2 

V S 1 V S 2 + V S 2 V S 3 + V S 3 V S 1 ~ V Sl V S a - V S2 V S 3 ~ % 3 %1 



n' 2 

_^L(A 2 -A 2 ). (3.9) 
8 

For a relatively large R\, A(Vo) vac < can happen; namely, the energy level of the sym- 
metric vacuum becomes lower than that of the broken one. Although such an EW vacuum 
could have a long lifetime that is larger than the age of the Universe and thus look viable, 
it is a highly nontrivial question whether the EW symmetry can be restored as the temper- 
ature increases. We thus exclude it from the investigation 1 . As mentioned in subsection 



1 In the viable MSSM baryogenesis scenario, the EW vacuum is metastable and sufficiently long-lived, 
and the charge-color-breaking vacuum is the global minimum. On the other hand, the energy level of the 
EW symmetric vacuum is higher than that of the broken one as in the usual scenario. In such a case, the 
successful EW symmetry restoration is possible H. 
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Figure 1: Left: A(V eS ) 
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(V r ff) vac as function of ra#±. Right: Av$ = v$ 



v s 



and AvSi — vs ; — v$ i7 i = 1,2,3, as a function of m H ±. We take tan/? = 1, A = 0.8, Ag = 0.3, 



v s = 500 GeV, v Sl 



vs 2 



1200 GeV, m 2 SSl 



4s 2 = ( 50 GeV ) 2 > "4s, = (200 GeV) 2 



3.1| , a large R\ is nothing but a large m H ±. Therefore, the maximally allowed value of 
rriff± can be derived from this vacuum condition. 

The signs of the second, fifth and last terms of Eq. |3.9| depend on the magnitudes 
of u's. In Fig. |], we plot various energy differences numerically. The left panel shows 
A(V c g) VSLC , where the one-loop corrections, V\ in Eq. (|3.1|) , have been included in the 
numerical calculation 2 . As an illustration, we take tan/3 = 1, A = 0.8, As = 0.3, vs = 500 
GeV, v.q, = vs. = = 1200 GeV, m% . = ml Q „ = (50 GeV) 2 , m| s = (200 GeV) 2 , 



VSi 

and A\ s 



SSi 



ss 2 



A\. We call this parameter set Case 1, which we will discuss in greater 
detail in Sec. [|. In this stable EW vacuum (A(T4ff) vac > 0) exists for 496 GeV 

< m H ± < 636 GeV. As mjj± decreases, the fifth term with a negative coefficient becomes 
dominant and then eventually results in A(V^,ff) vac < 0. On the other hand, as m H ± 
increases, A(V^f) vac < happens due to the contribution of the R\ term with a negative 
coefficient as indicated by Eq. (|3.9|) . 

In the right panel of Fig. [j], we plot Avs and Af^, i = 1, 2, 3, where Avs = vs — vs and 
Ai's, = — vs^ In the region where A(14ff) vac is small, \Avs\ and |AvsJ become large. 
As noticed in Ref . [^] , the locations of the symmetric and broken vacua at zero temperature 
may yield some information about the EWPT. Namely, the sizable |Aws| indicates that 
the singlet Higgs may be involved in realizing the non-MSSM-like EWPT. In such a case, 
since A(V r e fr) vac is small, the EW symmetry is expected to be restored at a relatively low 
temperature. 



4. Electroweak phase transition 

Here, we give the necessary ingredients for calculating the EWPT and describe the EWPT 



2 The behaviors of A(V e ff ) vac and A(Vo)vac are almost the same. 



-9- 



qualitatively. The one-loop effective potential at finite temperature takes the form 

Fi($ d ,$ u ,S,Si, 2j3 ;T) = I^^2 c aIb,f fe) , (4.1) 

with 

roc 

I BtF (a 2 ) = J dxx 2 In (l T e-^+^J , (4.2) 

where the subscripts of Ib,f{(i 2 ) denote boson (fermion). In order to reduce the com- 
putation time for the numerical integration in /^^(a 2 ), we will use the fitting functions 



employed in Rcf. [29] instead of Eq. ([4.2[). More explicitly, 

N 

/ B/ ( fl 2 ) = e- B ^ C ya", (4.3) 

71=0 

are used, where are determined by the least square method. For N = 40, \Ib,f{o- 2 ) — 
Ib,f(o. 2 )\ < 1CT 6 for any a, which suffices in our investigation. In the sMSSM, the structure 
of the tree-level potential is expected to be more important than the higher-order correc- 
tions unless a right-handed stop is lighter than top quark, which is the successful scenario 
of the MSSM EW baryogenesis. In the following, we exclusively explore a non-MSSM-like 
EWPT, i.e., the heavy stop case. Hence, we will not include the two-loop contributions, 
and neither will we perform the ring-improvements in the effective potential for simplicity. 

For the EW baryogenesis to work, the sphaleron process, which is active in the sym- 
metric phase, must be decoupled when the EWPT completes. In other words, the sphaleron 
rate in the broken phase should be less than the Hubble parameter at that moment. Con- 
ventionally, the sphaleron decoupling condition is cast into the form 

g > C, (4.4) 

where Tg is the temperature at which the EWPT ends, pe is defined by the vacuum 



expectation values of the two Higgs doublets at Te, i.e., pe = y pIC^e) + Pu(Te), and ( 
is an 0(1) parameter which depends on the profile of the sphaleron, etc. It is, however, 
a non-trivial task to evaluate Te explicitly since a full knowledge of bubble dynamics is 
required. We thus use the condition pc/Tc > C instead of Eq. (fO|), assuming that the 
supercooling is not too large, where Tq< is defined by the temperature at which the effective 
potential has two degenerate minima and p c is the VEV at Tq- To know the value of £ 
within 0(10%) accuracy, the sphaleron energy and zero-mode factors of the fluctuations 
around the sphaleron must be evaluated using the finite temperature effective potential. 
According to a recent study of the sphaleron decoupling condition with such an accuracy, 
£ ~ 1.4 in the MSSM [J29[j. In the NMSSM, the calculation of the sphaleron energy based 
on the tree-level potential has been done, and the sphaleron energy is found to be more or 
less the same as in the MSSM for a large portion of its parameter space Although 
a similar condition is expected for the sMSSM as well, we here adopt the rough criterion 
C = 1 for simplicity. 
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Figure 2: The patterns of the EWPT. 



4.1 Patterns of the EWPT 

As mentioned in Sec. ^ we will not consider the case of SCPV, which reduces the number 
of order parameters relevant to the EWPT to six. Now let us introduce 

P= (Pd,Pu,PS,PSi,PS2,PS 3 ) ( 4 -5) 

such that p(T = 0) = v = (v d ,v u ,vs,vs 1 ,vs 2 ,vs 3 ). 

As seen from Eq. ( |3.6| ), the effective potential in the ps direction has the form 

V(p s ) B -(R^+R^ps + ^mlpl + ^^p^ + ^Qs^QsA^Ps- (4-6) 

i 

Since the linear term in ps can exist in principle, the point (pd,p u ,ps) = (0,0,0) is not 
necessarily a local minimum. The coefficient of ps term can vanish only when pg x = ps 2 = 
3 . Therefore, we may categorize the EWPT into the following two types: 

Type A: (0,0, ps,ps 1 ,Ps 2 iPS 3 ) -> (Pd, Pu, PS, PSi, PS 2 , PS 3 ), 

Type B: (0, 0, 0, 0, 0, ) -4 (p d , p u ,ps, ps 1 , ps 2 , Ps 3 ) > 

where the barred quantities denote the corresponding VEV's in the symmetric phase. The 
phase transition pattern of the U(l)' symmetry may be even more diverse, and a detailed 
analysis of it is beyond the scope of the paper. In what follows, we will concentrate on 
the EWPT of Type A and B. We call it a Type AI EWPT if Av s > and a Type All 
EWPT if Avg < 0. These types of EWPT are pictorially shown in Fig. || 4 . In order to 
see the qualitative features of the EWPT, we consider a rather simplified case. In the case 
where the temperature is high compared to the mass of the particle in the loop, Ib,f((^) 



in Eq. (4.2) can be expanded in powers of a 2 = m 2 /T 2 as |3l[] 

I B (a 2 ) = -i- + ^a 2 - V) 3/2 " - ( In - - ~) + M,a« - • • • , (4.7) 

v ' 45 12 6 y J 32 V u B 2 J 384vr 2 v ; 

r , 2^ 77T 4 7T 2 , O 4 / O 2 3\ 7((3) fi , , 



360 24 32 I a F 2 J 384vr 2 



3 In the parameter space explored in this paper, sgn(7?i i 2) = +1 must be taken in order to be consistent 
with the positivity of the Higgs squared mass. 



4 



Of course there is additional axis of ps 3 which is not shown here. 
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where 



In ob = 21n4vr - 2^ E ~ 3.91, 



InaF = 21n7r - 2~/ E — 1.14, 



(4.9) 



and 7e(— 0.577) is the Euler constant and £(3)(~ 1.202) is a Riemann zeta function value. 
The relative errors of Ib,f{i 2 ) are less than 5% if a < 2.3 for bosons and a < 1.7 for 
fermions. 

In the multi-Higgs-doublet models, the (a 2 ) 3 / 2 term with a negative coefficient in 



Eq. (4.7) is crucial for the first order EWPT, as it gives rise to the potential barrier 
between the two degenerate minima. Such a cubic term originates from the zero frequency 
modes in the bosonic thermal loop. 

In the other limit, where the temperature is low compared to the mass, Ib,f{o. 2 ) can 
be expressed as [5pJ 



IbA* 2 ) - Ta 2 K 2 (a) ~ tJ\^ 2 ( 



1+ 8^ + 



(4.10) 



where K2(a) is the modified Bessel function. The relative errors of /^^(a 2 ) are less than 
5% if a > 2.6 for bosons and a > 2.3 for fermions. For the EWPT, the typical critical 
temperature is 0(100) GeV. Therefore, thermal effects of the squarks considered herein 
are exponentially suppressed, and do not play a significant role in realizing the first order 
EWPT. 



4.2 Type A 

We first begin with a simplified case. Suppose that g[ = and the VEVs of the secluded 
Higgs singlets are much larger than the others, which may correspond to the nMSSM-like 
limit. Although the EWPT in the nMSSM has already been studied in Refs. pSfl , we 
here give a brief review of it to get a qualitative feature of the EWPT in the sMSSM. Since 
the nMSSM does not possess the cubic self-interacting term of the singlet Higgs boson in 
the Higgs potential, the analysis of the EWPT becomes much simpler than that of any 
other singlet-extended MSSMs. 

Let us consider the EWPT in the subspace p = (pd,Pu,Ps) assuming pg i = vs v Here 
we consider the tree-level potential and the dominant temperature-dependent contribu- 
tions, which are proportional to T 2 . In the following, we show how the first order EWPT 
is realized without relying on the cubic term coming from the bosonic thermal loop as 
mentioned above. 

The effective potential here is reduced to 



V(p,T) = l -M\T)p 2 + l -m 2 p 2 



where 



(c + R x p 2 )ps+ ] ^p 2 p 2 s + ^p\ 



M 2 + GT 2 , 



M 2 (T) = m{ cos 2 (5 + m 2 , sin 2 /3 + QT = ,.,„ 
c = R\Vs 1 + R2VS2 > i?A = R\ sin (3 cos j3 
|A| 2 



A 



2 _ 92 + 9\ 



cos 2 2/3 + ^- sin 2 2(3, 



(4.11) 

(4.12) 
(4.13) 

(4.14) 
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and we have subtracted p-independent terms from V(p, T). The value of Q is given by the 
sum of the relevant couplings in the theory. For simplicity, the temperature dependence in 
the mixing angle (3 is neglected. In order to have a stable vacuum in the symmetric phase, 
m| > must hold. Using the tadpole condition of ps, ps can be written in terms of p 2 as 



Ps 



c + R\p 2 

|A|< 



mi + 



(4.15) 



which gives the trajectory of the minimum values of ps as a function of p. Plugging this 
back into the Higgs potential ( 4.11| ), we obtain 



V(p,T) = -M 2 (T)p 2 



(c + R x p 



2\2 



2(m 2 + \^p 



+ 



A 2 



-P 



(4.16) 



The problem is now reduced to a one dimensional EWPT analysis. The parameters Tc 
and pc for the first order EWPT are found to be 



T 2 



P 2 c 



F{pl) 



M 2 



F(v 2 



2 

w 



-m 2 s + 



2m| 



Ry 



\X\ 2 c 



2m| 



(4.17) 
(4.18) 



Here, F(p 2 ) is defined by 



F(p 2 ) = 2R X 



c + Rxp 2 
\M 2 



% ms+2~P J \ m s~^ 1 ~2~P' 




c + Rxp 2 
\M 2 



\ 2 p 2 . 



(4.19) 



From p c > 0, we find the condition for the first order EWPT |22| in terms of the model 
parameters: 5 



A < 



R\ 



\X\ 2 c 



2m 2 



(4.20) 



The physical implication of this condition becomes clearer if we expand the second term 



in Eq. ( 4.1 6| ) in powers of p 2 /m 2 s . To the sixth power of p, we obtain 



V(p,T) 



2mi 



+ I(m 2 (t)-^ 



+ - 



X 1 



Rx 




\Xfc 
Ami 



2m| 



Rx 



IX^c 
2m| 



(4.21) 



Condition flPoD requires that the coefficient of the /? 4 -term be negative. In other words, 
the roles of the negative cubic term and the positive quartic term in the usual scenario are 
now replaced by the negative quartic term and the positive sixth-power term. As already 
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pointed out in Ref . pjj , the form of the Higgs potential ( [4,21| ) is nothing but the SM Higgs 
potential with a dimension-six operator discussed in Refs. [j33|. Hence the nMSSM can be 
regarded as one UV completion of such an effective theory. 

It should be stressed that in order to make the first order EWPT stronger, m s should 
not be too large. This implies some constraints on the Higgs bosons whose masses come 
from m s . Another important implication is the following. From Eqs. ( [4.15 ) and (4.20), we 



find that to have a first order EWPT the difference between the singlet VEV in the broken 
phase and that in the symmetric phase must be larger than some critical value 



P 2 



,, 1 Xp 

\PS-Ps\ =J -2- [TT^f ' > , [Tj272-, (4.22) 

Trie 1 i A P In 1 1 _L A P 

s 1 + ^T y2m|l + i^- 

where ps = c/m s . Conversely, if p$ — ps, the singlet Higgs field will not play a significant 
role in realizing the first order EWPT, reducing to the MSSM-like EWPT. 

In the sMSSM case, because of the presence of the p^-term in the Higgs potential, the 
analytic formula for the first order EWPT is not as simple as the nMSSM case. However, 
it is the same mechanism at work. That is, the first order EWPT is possible due to the 
negative quartic term and the positive sixth-power term. 

Although Eq. ( |PcD can provide a good approximation to the nMSSM p3fl , it fails to do 
so quantitatively in the sMSSM due to the presence of the secluded singlet Higgs bosons. 
The precise value of ps strongly depends on pg i through the tadpole conditions in the 
symmetric phase. Indeed, it turns out that the above discussion is valid only qualitatively 
but not quantitatively. We will present the numerical results in Sec. ||. 

4.3 Type B 

As in the previous subsection, we focus only on the tree-level potential and the 0(T 2 ) cor- 



rections coming from the finite-temperature effective potential. As noted in subsection |4.1| , 
V can have an extremum at (p^, Pu, Ps, PSi > Ps 2 ) = (0> 0, 0, 0, 0) in Type B. It is thus useful 
to prameterize the />fields in terms of the five-dimensional polar coordinates 

Pd = z cos 5 cos 7 cos a cos /3, (4.23) 

p u = z cos 8 cos 7 cos a sin f3, (4.24) 

ps = z cos 5 cos 7 sin q, (4-25) 

ps 1 = z cos 5 sin 7, (4.26) 

ps 2 =z sin 5. (4.27) 

Using these variables, the effective potential at Tc takes the form 

V(z, T) = c 4 z 4 - c 3 z 3 + c 2 z 2 = c 4 z 2 (z - z c ) 2 , (4.28) 



where we have subtracted z-independent terms in Eq. (4.28), and 
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with 



C2 



222/ 22, 2 2 \ , 2272, 2 22, 2 2 

cgc^m^ + m 2Sl3 ) + m 5C(5 c 7 s Q + m Sl c s s^ + m S2 s s 



,2 2 



2 J2 2 „2 



,2 2 2 



,2 „2 



|A 5 | S 



-i?lC5S 7 C 7 S Q - R 2 SsCsC 1 S a - (i?i2 + -Ra s /5S 3 ) s <5 c <5 s 7 + 1 ~J J_ ( c i S 7 + ^PSa 



,2\ 2 



9? 



+-g c 2 5 c 2 c 2 a T 2 , 

C3 = -RAcfc^SaC^S^C^, 



C 4 



92 + 4 4 4 2 I^P 442/2 2 2, 2\, 
" c 5 c 7 C a C 2/3 ~7~ C S C 'y c a\ c a s P c P + S cJ + 



|A 5 | 2 



+ 



32 
9? 



2 2 2 



(4.30) 
(4.31) 

(4.32) 



where s a = sina,c a = cos a, etc and the angles (a, /3, 7, 5) are evaluated at Tc- Let us 
define c 2 = k% + k 2 T^. Then the critical temperature is given by 



A;2 V4c 4 



(4.33) 



In Type B, the magnitude of R\ is very important for the strong first order EWPT. 
This is the same as the usual scenario in which the first order EWPT is induced by the 
negative cubic term. Therefore, a large m#± is favored. We also note that a smaller C4 is 
preferred for a larger zq. This implies that a light Higgs boson is favored. We will quantify 
the statements here in Sec. ||. 

5. Numerical evaluations 



We give numerical results in this section. To this end, Tq and pc are determined using the 
one-loop effective potential at zero temperature, together with the fitting functions (^4.3^ of 



I b ,f{cl 2 ). First we discuss the CP-conserving case. The case of CP violation is argued in 



subsection |5.3| . Before showing the numerical results, we list the experimental constraints 
imposed in our numerical calculations. 

For a Higgs boson lighter than 114.4 GeV, we require 



i 2 < k(m Hi ) 



(5.1) 



where ^ = gn % zz / 'Th^zz an< ^ k 1S the 95% C.L. upper limit derived from the LEP exper- 
iments as a function of the Higgs boson mass. For the p parameter corrections, Ap < 
2.0 x 10~ 3 must be satisfied. As mentioned briefly below Eq. ( |3.1| ), the constraints of 
the Z' boson must be taken into account. The Z' boson mass (mz>) and the mixing be- 
tween Z and Z' bosons (azz') are constrained by the direct searches of the Z' boson and 
the EW precision measurements. The typical values are found to be mz> > 1000 GeV 
and azz> < 0(1O~ 3 ) p8fl . Since we here consider the azz> = case and at least one of 
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the VEV's of the secluded singlets is taken to be 0(1) TeV, the Z' constraints are easily 
satisfied. For the relevant SUSY particles, we impose m~± > 104 GeV and m^o > 46 GeV. 
In order to extract information about the doublet-singlet mixing effects, we define the 



MSSM fractions by [17 



?«, = (o<f>) 2 + (<>) 2 . «* = (<#>)', < 5 - 2 > 

where and are the orthogonal matrices which diagonalize the mass-squared 

matrices of the CP-even and CP-odd Higgs bosons, respectively. The parameter £ char- 
acterizes to what extent <j) (= H, A) comprises the doublet components. For £ = 1, <j) are 
purely composed of the doublets while £ = means that 4> is purely the singlet components. 

5.1 CP-conserving case 

We turn off CP violation in this and the next subsection. As an example we take 

Qn d = Qh u = 1, tan/3 = 1, |A S | = 0.3, A t = A b = 2m,- + // eff /tan/3, 

sgn(P As ) = sgn(P 1)2 ) = sgn(P 12 ) = +1, \m 2 SlSa \ = (200 GeV) 2 , 

m - q = m iR = m- bR = 1000 GeV, M 1 = 200 GeV, M 2 = M[ = 300 GeV, (5.3) 

where wig, and m^ R are soft SUSY breaking masses of the stop and sbottom, Mi, M 2 
and M{ are the gaugino masses associated with U(1)y, SU(2)l and U(l)', respectively. 
The quantity sgn(P.\) is fixed by 

sgn(P A ) = sgn ( m 2 H± - rr? w + — v 2 - Am 2 H± J , (5.4) 

where Am^ ± is the one-loop contribution to the charged Higgs boson mass. 

Since the mechanisms of the strong first order EWPT in Type A and Type B are dif- 
ferent from each other qualitatively, we will explore both parameter spaces. The following 
two cases are representative points for Types A and B, respectively. 

Case 1 : v s = 500 GeV, v Sl = v S2 = vs 3 = 1200 GeV, |A As | = \A X \, (5.5) 
Case 2 : v s = 500 GeV, v Sl = v S2 = 100 GeV, vs 3 = 1500 GeV, \A Xs \ = 1000 GeV. 

(5.6) 



As mentioned above, \A\\ is determined via the mass formula of mu±. The remaining 

l SSi I ' ' 



input parameters are A, mjy±, |m|c..|,i = 1,2, by varying which we search for the strong 



first order EWPT. 

In the left panel of Fig. ||, we plot pc, 2b, mjj 12 and as a function of m H ± 
in Case 1. As shown in Fig. |J the allowed region, where the prescribed EW vacuum is 
the global minimum, corresponds to the range, 496 GeV < m H ± < 636 GeV. It is found 
that there are two ranges in which pc/Tc > 1 is satisfied: 496 GeV < mjj± < 516 GeV 
(green-colored region) and 614 GeV < m H ± < 636 GeV (blue-colored region). The strong 
first order EWPT is possible only for the region with sizable |At>s|; |Au,g| > 110 GeV is 
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Figure 3: Dependence of pc, Tc, rriff 12 , rriA 1 and the MSSM fractions of H\, H 2 and A\ on m H ± 



in Case 1 with A 



0.8 and \m 2 SSi 



•ss 2 



(50 GeV)^ 



obtained in Case 1. In such a region, a low Tc is enough for the EW symmetry restoration 
as expected. We also plot the masses of Hi, Hi and A\. To be consistent with the strong 
first order EWPT, the mass of the lightest CP-even Higgs boson ra^ should be less than 
about 130 GeV. 

In the right panel of Fig. §, the MSSM fract ions of Hi, H2 and Ai are plotted. Devia- 
tions of £hi 2 from unity in the regions where the EWPT is strongly first order imply that 
the singlet Higgs bosons play a crucial role here. Since ^ 1 < 0.1, Ai is predominantly the 
singlet-like CP-odd Higgs boson. 

Next we examine the dependences of l^lsj and |m-|-5 2 | on pc/Tc- In Fig. |], we 

show pc, Tc, mn 1 and thai as a function of ^/l m isj> where \m 2 ss . | = Iw-lsJ = Iw-lsJ 
is assumed. As jw-I^J increases, pc/Tc decreases. Since a larger |m| s J gives a larger 
m 2 s via the tadpole condition, this tendency is understandable from the discussions given 
in subsection. [T^. Correspondingly, the Higgs boson masses which are mainly originated 
from m 2 s are constrained, leading to the mass bounds: < 135 GeV and m,A 1 < 143 

GeV to have the strong first order EWPT. 

The numerical results for Case 2 are shown in Fig. [B|. As done with Case 1, we 
plot pc, Tc, ™H]_;i an d (left panel) and £# 12 and (right panel) as a function of 
m#±. Similar to Case 1, due to the vacuum condition the allowed region is limited to 
the relatively small range 539 GeV < m H ± < 592 GeV. Within the interval, there are 
two regions where the strong first order EWPT is possible: 539 GeV < m H ± < 548 GeV 
(green-colored region) and 572 GeV < m H ± < 592 GeV (magenta-colored region). The 
former corresponds to Type AI EWPT and the latter to Type B EWPT. For Type B, 
as discussed in subsection ^L^, the magnitude of R\, and thus that of m H ± is crucial for 
the strength of the first order EWPT. As mjj± increases, pc/Tc becomes larger. From 
the right panel of Fig. ||, we find that the doublet-singlet mixing plays an essential role 
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Figure 4: Dependence of pc, Tc, mn 1 and m Al on \J\'m 2 ss \ 1 where IrafgJ = |w| s J = |wi| S2 | is 
assumed, in Case 1 with A = 0.8 and mjj± = 500 GeV. 
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Figure 5: Dependence of pc, 7b, 2 and (left panel) and 2 and (right panel) on 



m H ± in Case 2 with A = 0.8, \fng S | 



'ss 2 l 



(200 GeV) 2 



in realizing the strong first order EWPT as it should be. The mass limits of H\ and A\ 
consistent with pc/Tc > 1 are found to be, respectively, rriH 1 < 137 GeV and rriA 1 < 108 
GeV for this specific parameter set. 

5.2 Scan analysis 

In order to search for strong first order EWPT in the wider parameter space, we perform 
scans in the X-m H ± and \ \in 2 ss \-mH± planes, respectively. 
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I 2 I i 2 



We show the results for Case 1 in Fig. |(| In the left plot, we take \uKgg \ — i""s5 2 i 
(50 GeV) 2 and show the regions with sufficiently strong first order EWPT in the A-m#± 
plane. In the area surrounded by the solid black curves, the prescribed EW vacuum is 
the global minimum. The region below the magenta dashed line is excluded by the LEP 
95% C.L. exclusion limit, Eq. d57l|) . The green (blue) region corresponds to Type AI (All) 
EWPT. We also overlay the contours of the mass of the lightest Higgs boson, niH 1 = 140 
GeV (red), 160 GeV (green) and 180 GeV (orange). Around the A = 1 region, mH 1 = 140 
to 160 GeV is consistent with the strong first order EWPT. Here, we do not impose the 
perturbativity of A in the investigation. However, the maximally allowed mH\ satisfying 
Pc/Tc > 1 is expected to be larger than that of the MSSM. We will compare our results 



with the other models in subsection 5.4 



In the right plot, we take Case 1 with A = 0.8 in the w \m? ss , \-mn± plane and Im-sgJ = 

l 771 s , 5 2 l * s assumeci - The curves have the same meanings as the left plot except that the 
Higgs boson mass contours are drawn for A\. Here, we plot for mA 1 = 100 GeV (red), 
150 GeV (green), 200 GeV (orange) and 250 GeV (gray). It should be emphasized that 
rriA 1 < 250 GeV in order to be consistent with the strong first order EWPT. This feature 
differs from the MSSM. 

The numerical results of Case 2 is shown in Fig. ^. In the left panel, we take | = 

\m ss I = (200 GeV) 2 and show the regions with sufficiently strong first order EWPT 
in the A-m#± plane. The magenta region represents Type B EWPT. The meanings of 
the remaining curves are the same as in Case 1. The A dependence of Type B EWPT is 
expected from the analysis done in subsection L2. . Since C3 cx A and C4 oc A 2 in the potential 
04.28 ), pc/Tc becomes smaller as A increases. In both Type A and Type B EWPT, it is 
observed that mn 1 = 160 GeV supports the strong first order EWPT if A ~ 1. 
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In the right panel, we take A = 0.8 in the y\Trig S . \-mjj± plane, and IWggJ = |w|g 2 | = 
\m 2 ss . | is assumed. In Case 2, the dependence of | on the strength of the first order 
EWPT is smaller. This is because both vs 1 and vs 2 are taken to be small (100 GeV) so 
that the variation of m 2 s as a function of \rng S . | is much milder than that in Case 1. The 
m J 4 1 contours are drawn in the red curve (100 GeV), the green curve (150 GeV), and the 
orange curve (200 GeV), respectively. Asymptotically, tua 1 approaches around 214 GeV as 
|m| 5 .| increases. Therefore, there is no significant difference in the Higgs mass spectrum 
between Case 1 and Case 2. 

Here, we comment on the other cases. Although the magnitude of vs is relevant to the 
realization of the non-MSSM-like EWPT, so that it should be taken around 0(100) GeV, 
the maximal value of pc/Tc is not sensitive to its exact value. Similarly, the strength of the 
first order EWPT does not depend sensitively on the other unvaried parameters appearing 
in the tree-level Higgs potential. Therefore, the upper limits of the Higgs boson masses 
allowed for the strong first order EWPT would not change substantially. 

Even in the case of the non-MSSM-like EWPT, we could also argue that the light stop 
is lighter than the top quark. In Case 1 with the light stop, no significant enhancement on 
the strength of the first order EWPT is observed. It implies that the singlet Higgs bosons 
and the light stop do not contribute to the first order EWPT constructively. This may 
follow from the fact that the light stop gives a — p 3 -like term in the effective potential while 
the singlet Higgs bosons produces a — p 4 -like term. In Case 2 with the light stop, however, 
the stop contribution is constructive, leading to about 60% enhancement on pc/Tc- Since 
the right-handed soft SUSY breaking mass and the off-diagonal term in the stop mass 
matrix are taken to be small in the light stop scenario, the upper bound of mn\ obtained 
above is virtually unchanged. 



- 20 - 



5.3 CP- violating case 

We now discuss CP violation in this subsection. The relations between the CP-violating 



phases at the one- loop level are given by the one- loop tadpole conditions [19|: 
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(5.7) 
(5.8) 



where I g = lm.(XA q )/y/2, q = t,b. Here A q is the soft SUSY-breaking trilinear coupling, 
and / (m 2 ,m 2 ) is defined by 



/(mi, ml) 



mt — rrin 



(5.9) 



To see the CP-violating effect that the MSSM does not possess, we take It = h = so 
that I\ = 0. It follows from Eq. ( |5.&| ) that 
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-sm6 Sl s 2 , 



(5.10) 
(5.11) 



where 0sSi 2 = Arg(m| Si ). In the CP-violating case, the prefactors in Eqs. (5.10) and 



( 5.11 ) cannot be chosen arbitrarily. For instance, the following inequalities must hold for 
sin6» 5l 5 2 = 1: 



|m| lS2 |-J < |m| s J 



\m 



Si s 2 



< \m 

vs 



ss 2 \ 



(5.12) 



As we have discussed so far, small |^ss 12 | are favored for the strong first order EWPT 
in Case 1, which in turn imposes constraints on the magnitude of 5 | via Eq. ( 5.12j ). 
Since |m^ g \ appears in some of the diagonal elements of the neutral Higgs bosons, i.e., 
the (hs 12 )hs 12 ) an d (as 12 ,as 12 ) elements [19|, the Higgs boson masses mainly coming 



from those elements are lowered in accordance with \m s s |. 

In Fig. ||, we plot pc, Tq, and the three Higgs boson masses mn^i = 1,2,3, for Case 
1 with |A| = 0.8, m H ± = 630 GeV, \m 2 SSi \ = |mf<j 2 | = (50 GeV) 2 and |m| lS J = (20 GeV) 2 
(left panel); and for Case 2 with |A| = 0.8, m H ± = 590 GeV, and |m| s | = |m| S2 | = 
(200 GeV) 2 (right panel). In both cases, the CP-violating effect on the strength of the 
first EWPT is relatively mild. In particular, no significant dependence is observed in Case 
2. This is because the mass matrix elements of the CP-even and -odd mixing part are not 
sufficiently large to alter the Higgs boson masses significantly. 



Since |m| s | assumes a smaller value to be consistent with Eq. ( 5.12 ) in Case 1, some 
of the Higgs masses are lowered as seen in the left panel of Fig. [8|. In this case, Hi and H2 
are purely singlet-like Higgs bosons and may be challenging to detect at colliders. 
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Figure 8: Left: Case 1 with |A| = 0.8, m H ± = 630 GeV 
\m 2 SiS2 \ = (20 GeV) 2 . Right: Case 2 with |A 
(200 GeV) 2 . 
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Figure 9: Dependence of various quantities on the CP-violating phase #SiS 2 - Here we take A = 
0.8, m H ± = 675 GeV, v s = 560 GeV, v Sl = v Ss = 100 GeV, v S2 = 1000 GeV, \m 2 SSl \ = (334 GeV 2 ), 

KsJ = ( 106 Gev2 )> and \ m ks 2 \ = ( 20 ° Gey2 )- 



In contrast to Case 1, such a small \rrig g \ is not mandatory in Case 2 since vs 1 and 
vs 1 are small enough to satisfy Eq. ( p. 12 ). Therefore, there is no significant difference 
between the Higgs mass spectrum of the CP-conserving case and that of the CP-violating 
case. 

Depending on the choices of the Higgs VEVs, one can possibly find cases that have 
a significant dependence on 8s 1 s 2 - One example is shown in Fig. ||. Here we take |A| = 
0.8, m H ± = 675 GeV, v s = 560 GeV, v Sl = v Sa = 100 GeV, v S2 = 1000 GeV, \m 2 SSl \ = 
(334 GeV) 2 , |m| 52 | = (106 GeV) 2 , and |m| 1& | = (200 GeV) 2 . In this case, a stable EW 
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vacuum exists only for 0s 1 S2 ^ 90°) beyond which it shifts to a metastable one. It is found 
that pc/Tc < 1 at 0$ 1 s 2 = 0. As 0s 1 s 2 increases, however, pc/Tc gets enhanced and 
eventually pc/Tc > 1 is fulfilled for 0s 1 s 2 > 73.8°. This example clearly illustrates that 
CP violation in the model can enhance the first order EWPT. It should be noted that 
this behavior may not be observed in the MSSM, for the strong first order EWPT favors 
vanishing off-diagonal elements of the stop mass matrix, where the CP-violating phase 
resides. On the contrary, such a CP- violating phase can make the EWPT weaker in the 
MSSM §]. 

Finally, we briefly comment on the constraints from the electric dipole moment (EDM). 
Here, we only consider CP violation originating from the soft SUSY breaking masses m 2 SSi , 
m 5S , 2 an( ^ m li52 - ^ n sucn a case > the contributions to the EDM's of electron and neutron 



and so on are small enough to satisfy the current experimental limits [20]. 
5.4 Comparisons 

In this subsection, we compare our results with other models. The mass spectra having 
the strong first order EWPT are summarized in Table [l]. In the sMSSM, rriH 1 < 160 GeV 
and m J 4 1 < 250 GeV are required. However, there is no constraint on the stop mass. It is 
well-known that in the MSSM, in addition to the requirement of a light Higgs boson, the 
mass of the lighter stop must be smaller than that of top quark. According to Ref. 
rriH 1 < 127 GeV is required. An even severer bound raj < 120 GeV is found to be required 
by the strong first order EWPT. 

In the NMSSM, it is claimed that m,H 1 — 170 GeV is compatible with the strong first 



order EWPT [34]. In Ref 34], m,q = = 1000 GeV are taken, the two stop masses are 
assumed to be degenerate, and the perturbativity of A is taken into account. The sphaleron 
decoupling condition is pc/Tc > 1.3. Here, Tq is defined as the temperature at which the 
effective potential at the origin is destabilized in some direction. Such a temperature is 
always lower than the temperature that we define in this paper, rendering a larger pc /Te- 
la the nMSSM, the CP-odd Higgs boson must have an upper bound on its mass in 
order to realize the strong first order EWPT as discussed in subsection In Ref. p2| . 
such an upper limit is found to be < 250 GeV, which is approximately the same as 
the sMSSM. 

In the columns of NMSSM, nMSSM and UMSSM in Table 0, we leave a question mark 
for those cases where we are not aware of any literature that gives an upper bound on 
the Higgs boson mass consistent with the requirement of strong first order EWPT. As a 



reference, we cite the value of the CP-even Higgs boson in the nMSSM from Ref. [22]. It 
is found that the lightest CP-even Higgs boson with a mass of 130 GeV is consistent with 
the strong first order EWPT. The main differences between our work and Ref. p^ ] are 
the following: (1) The soft SUSY-breaking parameters are taken as rriq = mi = 500 GeV 
and At = 100 GeV. (2) The perturbativity of A (< 0.8) is imposed. (3) The sphaleron 
decoupling condition is pc/Tc > 1-3. When we naively adopt those three conditions in the 
sMSSM, the upper bound on itih 1 in the sMSSM approaches around 130 GeV. 



In the NMSSM |35| and the UMSSM |2J, it is found that the strength of the first 



order EWPT can be enhanced for a larger mass of the lightest Higgs boson. In contrast, 
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sMSSM 


MSSM 


NMSSM 


nMSSM 


UMSSM 


m Al 

mr 
tl 


< 160 GeV 

< 250 GeV 


< 127 GeV § 
< 120 GeV § < mt 


< 170 GeV || 

? 


? 

< 250 GeV H 


? 
? 



Table 1: The mass spectra for the strong first order EWPT in the various models. 



we do not observe such an effect in the sMSSM. 
6. Conclusions and Discussions 

We have investigated the possible parameter space where first order EWPT is possible and 
strong enough for successful EW baryogenesis in the sMSSM. We demonstrate two typical 
examples. In Case 1 where all of the secluded singlet Higgs VEV's are taken to be of 
the order of TeV, Type A EWPT is realized. This pattern of the EWPT is the same as 
in the nMSSM in the sense that the EWPT is possible by the non-cubic coupling with a 
negative coefficient. However, the dependences of the strength of the EWPT on m^± or 
A\ in both models are different from each other because of the U{1)' contributions and the 
secluded singlet sector. In Case 2 where two of the secluded singlet Higgs VEV's are taken 
to be 0(100) GeV, Type B EWPT is realized. In Type B, the parameters most relevant 
to the strong first order EWPT are m H ± and A. As m H ± increases, pc/Tc is enhanced. 
However, as A increases, the magnitude of pc/Tc reduces. The mechanism of strong first 
order EWPT in Type B is quite similar to the usual case where one has negative cubic and 
positive quartic terms. 

By scanning the parameters most relevant to the EWPT, we have obtained the typical 
Higgs mass spectrum that is consistent with the strong first order EWPT. In the sMSSM, 
it is found that mn x < 160 GeV and ttia 1 < 250 GeV must be satisfied. Unlike the MSSM, 
the light stop mass is not necessarily smaller than the top quark mass. 

We have also worked out the impact of the CP- violating phase on the strength of 
the first order EWPT. In most of the parameter space, the dependence of pc/Tc on 
6siS 2 is quite mild. However, according to the tadpole conditions for the CP-odd Higgs 
bosons, \mg S |, |m<!s 2 | and |w| iS2 |, vs and vs ia cannot be freely chosen. This leads to the 
constraints on the Higgs mass spectrum, especially for Case 1. 

Our numerical study suggests that to have non-MSSM-like EWPT, the singlet Higgs 
VEV's, particularly vs, in the broken phase and the symmetric phase must be significantly 
different from each other. For the typical parameter sets, \Avs\ > 100 GeV must be 
satisfied. This condition is almost temperature-independent in our analysis. In principle, 
Avs can be derived provided that the soft SUSY-breaking masses are known, or more 
precisely, once the global structure of the Higgs potential is completely determined. The 
determination of a sizable |A-us| at zero temperature from collider experiments may be 
evidence of strong first order EWPT in the singlet-extended MSSM. 
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